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In this paper, we give a positive answer to the open question: Can there exist 4 limit cycles in 
quadratic near-integrable polynomial systems? It is shown that when a quadratic integrable 
system has two centers and is perturbed by quadratic polynomials, it can generate at least 
4 limit cycles with (3, 1) distribution. The method of Melnikov function is used. 
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1 Introduction 

The well-known Hilbert's 16th problem is remained unsolved since Hilbert [1] proposed the 
23 mathematical problems at the Second International Congress of Mathematics in 1990. Re- 
cently, a modern version of the second part of the 16th problem was formulated by Smale [2] , 
chosen as one of the 18 challenging mathematical problems for the 21st century. To be more 
specific, consider the following planar system: 

^ = Pn(x, y), ^ = Q n (x, y), (1.1) 

where P n (x,y) and Q n (x,y) represent n th -degree polynomials of x and y. The second 
part of Hilbert's 16th problem is to find the upper bound H(n) < n q on the number of 
limit cycles that the system can have, where q is a universal constant, and H(n) is called 
Hilbert number. In early 90's of the last century, Ilyashenko [3j and Ecalle [4j proved the 
finiteness theorem pioneered by Dulac, for given planar polynomial vector fields. In general 
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the finiteness problem has not been solved even for quadratic systems. A recent survey 
article [5] (and more references therein) has comprehensively discussed this problem and 
reported the recent progress. 

If the problem is restricted to the neighborhood of isolated fixed points, then the question 
is reduced to studying degenerate Hopf bifurcations, which give rise to fine focus points. In 
the past six decades, many researchers have considered the local problem and obtained many 
results (e.g., see [6-12]). In the last 20 years, much progress on finite cyclicity near a fine 
focus point or a homoclinic loop has been achieved. Roughly speaking, the so-called finite 
cyclicity means that at most a finite number of limit cycles can exist in some neighborhood 
of focus points or homoclinic loop under small perturbations on the system's parameters. 

In this paper, we particularly consider bifurcation of limit cycles in quadratic systems. 
Early results can be found in a survey article by Ye [13] . Some recent progress has been 
reported in a number of papers (e.g., see [HI US]). For general quadratic system (11. ip 
(n — 2), in 1952, Bautin [5] proved that there exist 3 small limit cycles around a fine focus 
point or a center. After 30 years, until the end of 1970's, concrete examples were given 
to show that general quadratic systems can have 4 limit cycles [16j [17], around two foci 
with (3, 1) configuration. Since then, many researchers have paid attention to integrable 
quadratic systems, and a number of results have been obtained. A question was naturally 
raised: Can near-integrable quadratic systems have 4 limit cycles? A quadratic system is 
called near-integrable if it is a perturbation of a quadratic integrable system by quadratic 
polynomials. On one hand, it is reasonable to believe that the answer should be positive 
since general quadratic systems have at least 4 limit cycles; while on the other hand, near- 
integrable quadratic systems have limitations on their system parameters and thus it is more 
difficulty to find 4 limit cycles in such systems. In fact, this is still an open problem after 
another 30 years since the finding of 4 limit cycles in general quadratic systems. 

The study of bifurcation of limit cycles for near-integrable systems is related to the so 
called weak Hilbert's 16th problem [18J, which is transformed to finding the maximal number 
of isolated zeros of the Abelian integral or Melnikov function: 

M(h,5)=i Q n dx- P n dy,lb2 (1.2) 

JH(x,y)=h 

where H(x,y), P n and Q n are all real polynomials of x and y with degH = n + 1, and 
max{degP n , degQ n } < n. The weak Hilbert's 16th problem is a very important problem, 
closely related to the maximal number of limit cycles of the following near-Hamiltonian 
system [T9] : 

dx dH(x,y) dy dH(x,y) . 

~dt = — dy — + £ P^ x ^y^ ~fa = fa — + eq n {x, y),lb3 (1.3) 

where H(x, y), p n (x, y) and q n (x, y) are polynomials of x and y, and < e <C 1 is a small 
perturbation. 

General quadratic systems with one center have been classified by Zoladek [20J using a 
complex analysis on the condition of the center, as four systems: Q^ v - the Lotka-Volterra 
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system; - Hamiltonian system; Q% - reversible system; and Q4 - codimension-4 system. 
In 1994, Horozov and Iliev [21] proved that in quadratic perturbation of generic quadratic 
Hamiltonian vector fields with one center and three saddle points there can appear at most 
two limit cycles, and this bound is exact. Later, Gavrilov [22] extended Horozov and Kiev's 
method to give a fairly complete analysis on quadratic Hamiltonian systems with quadratic 
perturbations. Quadratic Hamiltonian systems, with at most four singularities, can be clas- 
sified as three cases [22J : (i) one center and three saddle points; (ii) one center and one saddle 
point; and (iii) two centers and two saddle points. In [22], Gavrilov showed that like case 
(i), cases (ii) and (iii) can also have at most two limit cycles. Therefore, generic quadratic 
Hamiltonian systems with quadratic perturbations can have maximal two limit cycles, and 
this case has been completely solved. 

For the Qf reversible system, there have been many results published. For example, 
Dumortier el al. [23] studied a case of Q§ system with two centers and two unbounded 
heteroclinic loops, and presented a complete analysis of quadratic 3-parameter unfolding. It 
was proved that 3 is the maximal number of limit cycles surrounding a single focus, and 
only the (1, 1) -configuration can occur in case of simultaneous nests of limit cycles. That is, 
3 is the maximal number of limit cycles for the system they studied [23] . Later, Peng [23] 
considered a similar case with a homoclinic loop and showed that 2 is the maximal number 
of limit cycles which can bifurcate from the system. Around the same time, Yu and Li [25] 
investigated a similar case as Peng considered but with a varied parameter in a certain 
interval, and obtained the same conclusion as Peng's. Later, Iliev et al. |26j re-investigated 
the same case but for the varied parameter in a different interval (which yields two centers) 
and got the same conclusion as that of [23], i.e., 3 is the maximal number of limit cycles 
which can be obtained from this case. Recently, Li and Llibre [2Z] considered a different case 
of Qf system which can exhibit the configurations of limit cycles: (0,0), (1,0), (1,1) and 
(1,2). Again, no 4 limit cycles were found. In order to explain why the above authors did 
not find 4 limit cycles from the Qf reversible system, consider the system with quadratic 
perturbations, which can be described by [23] 

x = -y + ax 2 + by 2 + e^x + /j 2 xy), 

IbA (1.4) 

y = x (1 + cy) + 5^3 x, 

where a, b, c are real parameters, /Uj, i — 1, 2, 3 are real perturbation parameters, and 
< £ < 1. When 5 = 0, system (??) e= o is a reversible integrable system. It has been noted 
that in all the cases considered in [231 122 123 121] , the parameters a and c were chosen as 
a = - 3, c = - 2, but with 6 = 1 in [23]; b= -1 in [23], b e (- 00, -1) U (-1, 0) in [25], 
and b G (0, 2) in [26]. In these papers, complete analysis on the perturbation parameters 
was carried out with the aid of Poincare transformation and the Picard-Fuchs equation, but 
it needed to fix all (or most of) the parameters a, b and c. This way it may miss opportunity 
to find more limit cycles, such as possible existence of 4 limit cycles. As a matter of fact, 
for the cases considered in [251 EE] , a simple scaling on the parameter b {b 7^ 0) can be used 
to eliminate b. So, suppose the non-perturbed system (??) e=0 has two free parameters and 
let us consider the 2-dimensional parameter plane. Then, all the cases studied in the above 
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mentioned articles are special cases, represented by just a point or a line segment in the 
2-dimensional parameter plane (see more details in Section 2). It has been noted that a 
different method was used in [27] with Melnikov function up to second order, but no more 
limit cycles were found. 

It should be mentioned that Zhang [28] has proved that the possible cycle distributions 
in general quadratic systems with two foci must be (0, redistribution or (1, ^-distribution, 

1 = 0, 1, 2, 3, • • •. So far, no results have been obtained for i > 4. This result also rules 
out the possibility of (2, 2)-distribution. It is conjectured that at most 3 limit cycles can 
exist around one focus point. The problem of bifurcation of 3 limit cycles near an isolated 
homoclinic loop is still open. 

In this paper, we turn to a different angle to consider bifurcation of limit cycles in 
quadratic near-integrable systems with two centers. We shall leave more free parameters in 
the integrable systems, so that we will have more chances to find more limit cycles. The basic 
idea is as follows: we first consider bifurcation of multiple limit cycles from Hopf singularity, 
which does not need to fix any parameters, and use expansion of Melnikov function near 
centers to get such limit cycles as many as possible. This leads to determination of a maximal 
number of parameters. Then, for the remaining undetermined parameters, we compute the 
global Melnikov function to look for possible large limit cycles. Indeed, although, due to 
the complex integrating factor in the analysis, we are not able to give a complete analysis 
for classifying the perturbation unfolding, we do get a positive answer to the open question 
of existence of 4 limit cycles in quadratic near-integrable systems. In particular, we will 
show that perturbing a reversible, integrable quadratic system with two centers can have at 
least 4 limit cycles, with (3,1) distribution, bifurcating from the two centers under quadratic 
perturbations. 

The rest of paper is organized as follows. In Section 2, we give a different classification 
in real domain for quadratic systems with one center, and compare it with that given by 
Zoladek [20J. Also, we use our classification to present a simple summary on some of the 
existing results for the reversible near-integrable system. Section 3 is devoted to the analysis 
on bifurcation of small limit cycles from Hopf singularity. In Section 4, we show how to find 
large limit cycles bifurcating from closed orbits to obtain a total of 4 limit cycles. Finally, 
conclusion is drawn in Section 5. 

2 Classification of generic quadratic systems with at 
least one center 

In this section, we give a different classification in real domain for quadratic systems with 
a center, which is consistent with the Hamiltonian systems considered in [2T| [22] . We start 
from the following general quadratic system: 
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d,Z\ 2 2 

~TT = c 100 + c 110 z l + c 101 z 2 + C120 + Clll ^1 z 2 + C102 ^ 2 , 

f 165 (2.1) 

2 _ i i , 2 , ,2 

— TT — C200 + c 210 ^1 + C201 z 2 + C22O ^1 + C211 2i 2^2 + C202 

at 

where cafe's are real constant parameters. It is easy to show that this system has at most four 
singularities, or more precisely, it can have 0, 2 or 4 singularities in real domain. In order for 
system (??) to have limit cycles, the system must have some singularity. In this paper, we 
assume that system (??) has at least two singularities. Without loss of generality, we may 
assume that one singular point is located at the origin (0,0), which implies cioo = C200 — 0, 
and the other at (p,q) (p 2 + q 2 7^ 0). Further assume the origin is a linear center. Then 
introducing a series of linear transformations, parameter rescaling and time rescaling to 
system (??) yields the following general quadratic system: 

dx 2 

— = y + a x xy + a 2 y , 

f M (2.2) 

u y 1 2 1 1 2 

— = - x + x + a 3 xy + a A y , 
at 

which has a linear center at the origin (0,0) and another singularity at (1,0). 

In order to have the origin of system (??) being a center, we may calculate the focus 
values of system (??) and find four cases under which (0, 0) is a center, listed in the following 
theorem (here we use Zoladek's notation in our classification). 

Theorem 1.1 The origin of (11) is a center if and only if one of the following conditions 
is satisfied: 

-Reversible system: a 3 = a 2 = 0, under which system (11) becomes 

dx 

— = y + a t xy, 

f im (2.3) 

d v , 2 , 2 

— = - x + x +a 4 y , 
at 

with 

{center if a± < — 1 , 

saddle point if a± > — 1. 

- Hamiltonian system: a 3 = a\ + 2 a 4 = 0, under which system (11) is reduced to 

^ = y + a 1 xy + a 2 y 2 , 

% 1 2 '» (") 

Tt = ~ x + X ~2 aiV - 

with 

( center if a-i < — 1 , 

T, 0) being a { 

saddle point if a± > — 1. 
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Qg - Lokta-Volterra system: a 2 = 1 + a 4 = 0, under which system (11) becomes 

dx 

— = y + atxy, 

f IW (2.5) 

— = -x + x 2 + a 3 xy-y 2 , 
dt 



with 



focus if (ii < — (1 + j a 2 



4 "3 )i 

1 „2 



'1, 0) being a ^ node if — (1 + | a| ) < ai < — 1, 

saddle point if ai > — 1. 

Q4 - Codimension-4 system: 

a 3 - 5 a 2 = ai - (5 + 3 a 4 ) = a 4 + 2(1 + a^) = 0, /WO (2.6) 

under which system (11 ) can be rewritten as 

= V ~ (1 + 6 aj) x y + a 2 y 2 , 
dt " Ml (2.7) 

— = — x + x 2 + 5a2xy — 2(1 + a\) y 2 , 

(i/L 

with (1,0) being a node for a 2 7^ 0. 

Remark 1.2. There is one more case found from the above process, defined by the following 
conditions: 

a 3 - 5 a 2 = ai - (5 + 3 a 4 ) = 3 (a 4 + 2) (a 4 + l) 2 - (5 a 4 + 6) a 2 = O.lbU (2.8) 

We will show later in this section, when we compare our above real classification with the 
complex classification given by Zoladek [20] . that the case defined by (??) actually belongs 
to the Qf -reversible system. 

Proof. Necessity is easy to be verified by computing the focus values of system (??) associ- 
ated with the origin. Some focus values will not equal zero if the condition is not satisfied. 

For sufficiency, we find an integrating factor for each case when the condition holds. 
For the Q% - Hamiltonian system (??), we know that the integrating factor is 1, and the 
Hamiltonian is given by 

H(x, y) = X - [x 2 + y 2 ) - ± x 3 + 1 a x x y 2 + 1 a 2 y 3 , IbU (2.9) 

which is exactly the same as that given in [211 E2] . 

For the Q% - reversible system (??), the integrating factor is 

7 = |1 + a 1 x\~^T 1 , IbU (2.10) 
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and the first integral of the system is given by 



F(x, y) — — sign(l + a^x) |1 + a\x\ a i 



2 ^ (1 + ai — a 4 ) (1 + 2 a 4 x) 



x 



.Iblh 
(2.11) 



a 4 (ai — a 4 ) (ai — 2 a 4 ) a 4 — a 4 
For the Q% v - Lokta-Volterra system (??), we find the integrating factor to be 
7 = \g{x^y)\~ x , where g(x, y) = (1 + a x x) [(x- l) 2 + a 3 (x - 1) y - (1 + a 1 ) y 2 \ , Z616 (2.12) 
and the first integral of the system is 

- sl S n ^( x ^y)) | 2 i n |i +aia; | +ai i n (l+ai)2/ 2 — a 3 2/(^— 1) — (^— 1) 2 | 
2 ai(l + ai) L 



F(x,y) 



+ 



2 ai ci3 (x— 1) 



^[a2+4(l+a 1 )] (x-l)2 

when a 2 + 4 (1 + a%) > 



tanh 1 



q 3 (a;-l)-2(l+ai)y 
v /[a§+4(l+a 1 )] (*-!)= 



} 



sign(l + aix) 
2ai(l + ai) 



j^2 In |l+aix| +a 4 In (l + ai)y 2 — a 3 y(x— 1) — (x— 1)' 



Z617 



2 ai 03 (ir— 1) 



y/[-a*-4(l+ai)] (x-1) 

when a 2 + 4 (1 + ai) < 0. 



tan 



-i 



a 3 (a:-l)-2(l+ai)y 
v /[_ a 2_ 4 ( i +ai) ] ^Zi)j 



(2.13) 



Finally, for the Q 4 - codimension-4 system (??), we have 

1 = \ g (x,y)\' 5/2 , where g(x, y) = 1-2 (1+2 a 2 ) x-2 a 2 y+(l+4 a\) {x+a 2 y) 2 , lbl8 (2.14) 

and the first integral of the system is equal to 

sign(#(x,y)) 



F(x,y) 



12 a 6 



\g(x,y)r 3 / 2 f(x iy ),lbl9 



where 



f(x, y) = — (1 + a\) + 3 (x + a 2 y + 2 a\ x) 1 + a 2 , — (1 + 3 a 2 ,) (x + a 2 2/) 
+ (1 + 3 a 2 ) (1 + 4 a 2 ) (x + a 2 y) 3 . 



(2.15) 



/620 (2.16) 



The proof is complete. □ 

Note that among the four classifications of the integrable system (??), the first three 
classified systems (??), (??) and (??) have two free parameters, while the last system (??) 
has only one free parameter. 

Remark 1.3. We now show that our classification in Theorem 1.1 is equivalent to that given 
by Zoladek [20J. The general quadratic system considered in [20] is given in the complex 
form: 

dz 

= (i + \) z + Az 2 + Bzz + Cz 2 ,lb21 (2.17) 

at 
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A = 


B = 0, 






Qf '■ 


A = 


2A + B 


= 0, 




Of: 


A = 


lm(AB) 


= lm(B 3 C) 


= lm{A 3 C) = 0, 


Qa- 


A = 


A-2B 


= \C\-\B\ 


= 0. 



where z = x + i y, and A, B and C are complex coefficients. It has been shown in [20J that 
the point z = is a center if and only if one of the following conditions is fulfilled: 



W22 (2.18) 



In the following, we first use real differential equation to give a brief proof (different 
from Zoladek's [20J ) , and then show that our classification is equivalent to Zoladek's when 
system (??) is assumed to have a non-zero singularity. To prove this, let 

A = A 1 + iA 2 , B = B 1 + iB 2 , C = d+zC 2 , {i 1 = - 1), 

and then rewrite the complex equation (??) in the real form: 

dx 



t \x + y + (A 1 + B l + C 1 )x 2 + 2 (A 2 -C 2 )xy- (A 1 -B 1 + C x ) y 2 , 
dt W23 (2.19) 

- x + Xy - (A 2 + B 2 + C 2 ) x 2 + 2 {A x -C 1 )xy + (A 2 - B 2 + C 2 ) y 2 , 



dt 

where y —> —y has been used. Letting A = yields the focus value v = 0. Then, it is 
easy to find the first focus value (or the first Lyapunov constant) as 

V! = - A x B 2 - B t A 2 = - lm(AB).lb24 (2.20) 

Letting v± = results in lm(AB) = 0, which gives 

B A 

B 2 = — , under the assumption of A\ ^ 0.lb25 (2-21) 

A\ 

(The degenerate case A\ = can be similarly analyzed and the details are omitted here.) 
Then, we apply our Maple program (e.g., see [29]) to system (??), with the conditions A = 
and (??), to obtain 

-f(A x -2B x ) -ff 3 -fU ~fh 



V 2 = TTTi > V 3 = TS , V A = — — A7 , V 5 



3A 3 ' 6 216Af 4 9720AJ' 466560A?' 
where 

f = Bi {2A 1 + B l ){C 2 A 3 1 + 3C 1 A 2 A\ - ?>A\C 2 A X - C X A\), 
and fs, /4, etc. are polynomials of A x , A 2 , C x , C 2 and B x . Letting / = 0, i.e., 

B x = or 2 A x + B x = or C 2 A\ + 3 C x A 2 A\-3 A 2 2 C 2 A x — C X A\ — lm(A 3 C) = 

yields v 2 — v 3 = • • • = 0. 
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Indeed, By = implies B 2 = due to the condition (??), and so B — 0. Thus, we 
obtain A = B = 0, corresponding to the Q% v case. 

For the condition 2 A 1 +B l = 0, it follows from (??) that 2 A 2 -B 2 = 0, i.e., 2 A+S = 0, 
which plus the condition A = gives the Q% case. 

The third condition \m(A 3 C) = 0, with A = and Im(AB) = 0, corresponds to the 
Qf case. Further, it is easy to show that under the condition lm(AB) = 0, lm(A 3 C) = 
and Im(£> 3 C) = are equivalent. Thus, the conditions A = lm(AB) = lm(B 3 C) = also 
applicable for this case. So for this case, either lm(A 3 C) = or lm(B 3 C) = is needed, 
but not both of them. In the following, we show one more case to join this case, leading to 
both the two conditions being needed. 

Note that there is one more condition A 1 = 2 Bi which renders v 2 = 0. Letting 
Ai — 2 B-y, and so A 2 = — 2 B 2 (see (??)), implying that A — 2 B = 0. Under the condition 
A = 2 B, v\ = v 2 = 0, and the other focus values become 

25 



V3 

V5 = 



6480 



2 s ) 

^ = (Cl + C\ - B\ - Bl){C 2 B\ - ZC X B\B* - ZC*B X B\ + dS 2 3 ), 

o 

\hB\ + 5855^ + 60( J B 1 C 1 + B 2 C 2 ) - 196(C 1 2 + C 2 2 )] , 

648(75f + 124S 2 S 2 2 + 1557S 2 1 )-3(961fi 2 C 2 -76805i5 2 C , iC' 2 + 202345 J B 2 C 2 ) 

+ 576fiiCi(106fi 2 + 307S 2 ) + 288B 2 C 2 (371S 2 + 773fi 2 ) 
- 3(4801 J B 1 2 C 2 2 + 206185 J B 2 C 1 2 ) - 80688(C1 2 + C2 2 )(B 1 C 1 + B 2 C 2 ) 
+ 86144(C1 2 + C2 2 ) 2n 



Hence, under the conditions X = A — 25 = 0, there are two possibilities such that v% = 
t> 4 = • • • = 0. The first possibility is 

d + Cl-Bl-Bl i.e., |C|-|fl|=0, 

which is one of the conditions given for the case (see (??)). 
The second possibility is given by the condition: 

C 2 B\ -3d B\ B 2 -3C 2 B x B\ + C x B\ = lm(B 3 C) = | lm(A 3 C) = 0, lb2Q (2.22) 

due to A = 2 B. Since these conditions can be included in the conditions A = Im(AB) = 
lm(B 3 C) = lm(A 3 C) = 0, this possibility belongs to the Q§ case. 

The remaining task is to show that the conditions classified in (??) are sufficient. This 
can be done by finding an integrating factor for each case. For brevity, we only list these 
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7 = < 



integrating factors below (while the lengthy expressions of the first integrals are omitted): 

1 + 4 (A 2 x - A l y) + 4 (Ad + A 2 d - 2A 1 A 2 ) x y 
+ [(A, + d) (A 1 - 3d) + (A 2 + C 2 ) (5A 2 - 3d)] x 2 
+ [(A 2 + C 2 ) (A 2 - 3C7 2 ) + (A + d) (5A X - 3d)] y 2 
+ 2 (A 2 + A 2 - C 2 - C 2 ) [(A 2 + C 2 ) x 3 - (A, + d) 
-(Ai-3d)z 2 y + CA 2 -3d):ry 2 ] " J 



1 - 2 {A 1 - d) 2/| 
1 - 4(S 2 x + B lV ) + 2{B\ + B 2 ) {x 2 + y 2 ) 



+ 2{B l C l + B 2 C 2 ) (:r 5 



i(B 1 C 2 -B 2 C 1 )xy 



-5/2 



for Q3 , 

for Qf, 

for Qf, 

for Q4. 



Ib27 (2.23) 



For the integrating factors of degenerate cases (e.g., A\ — C\ = 0), one can easily find them. 

Next, compare the classification listed in (??) with ours given in Theorem 1.1. First, 
consider the Q% v case. Letting A = B 1 = B 2 = in (??) yields 



(IT 

— = y + (A 1 + C 1 )x 2 + 2 (A 2 -C 2 )xy- (A, + d) y\ 
^L = - x -(A 2 + C 2 )x 2 + 2(A l -C 1 )xy+(A 2 + C 2 )y 2 . 



W28 



(2.24) 



Then, let 



k = tan(#), and so sin(6') 



k 



:, C0s(#) 



1 



=, Z629 



VrTk 2 ' VTT¥ 

where k is solved from the following cubic polynomial: 

P 1 (k) = (A 2 + C 2 ) k 3 + (Ai - 3 d) k 2 + (A 2 - 3 d) k + A 1 + C 1 = 0./630 



(2.25) 



(2.26) 



This cubic polynomial at least has one real solution for k, which gives the slope of the line 
on which a second fixed point is located, k = if A\ + d = 0, otherwise, k ^ 0. Let k be 
a real root of Pi(k), i.e., -Pi (A;) = 0. 

Further, introducing the linear transformation (rotation): 



x = cos(#) u — sin(6>) v, y = sin(#) u + cos(#) v, Z631 



into (?? yields 



^ = y + mi20 + mm x y + mio2 y 2 , 
at 

dy o 2 

— = - x + m 220 x + m 211 xy + m 202 y , 
at 



lb32 



(2.27) 



(2.28) 
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where 

m 120 = -m W2 = (l + kY^P^k) = 0, 

m 220 = - m 202 = (1 + k 2 )- 3/2 [(A + d) t - (A 2 - 3 C 2 ) fc 2 + - 3 d) - A 2 - d 
mm = -2(1+ 1 2 )" 3/2 [(^i " Ci) I 3 - (A 2 + 3 d) I 2 + + 3 d) fc - A 2 + d 
man = 2(1+ F^ 3 / 2 [(A 2 - d) t + (A 1 + 3 d) £ 2 + (A 2 + 3 d) fc + A 1 - d 

Suppose m 220 7^ 0. Then, introducing x = m 220 x, y = m 220 y into (??) results in 

dx m m _ 

-77 = y + zy, 

f m220 Z633 (2.29) 

^ - , -2 , m 2ii _ 2 

at m 2 2o 

which is identical to (??) as long as letting 01 = and a?. = 22211. This shows that 

v I & & 1 m 2 20 ° »™220 

the four parameters A 1: A 2 , d and C 2 are not independent. Thus, alternatively, we may 
simply take k = (which renders the second singularity of (??) on the rr-axis), yielding 
d = ~A X . Thus, (??) becomes 

^ = y + 2(A 2 -C 2 )xy, 

^L = - x -(A 2 + C 2 )x 2 + 4A 1 xy + (A 2 + C 2 )y 2 . 
at 

Suppose A 2 + C 2 7^ 0. Introducing x = — (A 2 + C 2 ) x, y = — (A 2 + C 2 ) y into the above 
equations we obtain 

^ _ - _ 2(A 2 -C 2 ) 

f_ /&34 (2.30) 

ct y 2 A A. 1 2 

7^ = _x + x ~~ A2+C2 xy ~~ y ' 

which is identical to (??) if letting ai = ~ 2 4 ^~^' 2 - > and a 3 = +^f^ ■ In the following, we 
will use this simple approach for other cases. 

For the Q% case, substituting A = 0, Bi = — 2 A 1 and B 2 = 2 A 2 into system (??) 
results in 

— = y -(A 1 -C 1 )x 2 + 2(A 2 -C 2 )xy-(3A 1 + C 1 )y 2 , 
^ = -x-(3A 2 + C 2 )x 2 + 2(A 1 -C\)xy- (A 2 - C 2 ) y 2 . 

Further, taking d = A\ in the above equations gives another singularity on the x-axis, and 
introducing x = — (3 A 2 + C 2 ) x, y = — (3 A 2 + C 2 ) y into the resulting equations yields 



P. Yu & M. Han 



4 limit cycles in near- quadratic systems 



12 



which is identical to (??) if we set ai = - |^! [ _~^P 2 ' ) and a 2 = 3 ^ A ^ C ■ 

For the Qf reversible case, it follows from [20] that all the coefficients A, B and C are 
real, and thus we obtain the following real form from the complex system (??) 

(IjX r) 7 9 

_ = -y + ax + by , 

f lb36 (2.32) 

_ = x + cxy, 
at 

where 

a = A 1 + B 1 + C 1 , b = B 1 -A 1 -C 1 , c = 2A 1 -2C 1 . 
Suppose b 7^ 0. Then, introducing x = by, y = bx into (??) results in 

dx _ c 

— = y+- X y, 

$. b mi 2.33 

dy _ A _ 2 . a _ 2 v ' 

Tt = ~ x + X + b y > 

which is identical to (??) if 

c 2{A l -C l ) J a A 1 + B 1 + C 1 

ai = b = B x - Al - Cl ^ ai= b = B.-M-C,' 

For the last Q 4 case, under the condition X = A — 25 = 0, by setting C\ = — 3 B\ 
(which renders a non-zero singularity on the x-axis) in (??) we obtain 

dx 

— = y-2(2B 2 + C 2 )xy + 2B iy 2 , 
dt 

^ = - x + (B 2 - C 2 ) x 2 + 10 B l x y - (3 B 2 - C 2 ) y 2 , 
dt 

Suppose B 2 — C 2 7^ 0. Then, introducing x = (B 2 — C 2 ) x, y = (B 2 — C 2 ) y into the above 
equations yields 



d % 2{2 B 2 +C 2 ) i 2B 1 -2 



—r = V p n x y + -5 — y 

* B *- C2 "'- C2 (638 (2.34) 

Comparing the coefficients of the above system (??) with our system (??) results in 

which in turn implies that 03 — 5 a 2 = a\ — (5 + 3 04) = 0, and 

a 4 + 2 (1 + a 2 2 ) = 8B { lt%? = C ' + { S-c!r Bl = °» for |C| - = 0. 
The above conditions are the exact conditions given in (??) for the Q4 case. 
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Finally, we turn to the conditions given in (??). It follows from (??) that 

3 (04 + 2) (04 + l) 2 - (5 a 4 + 6) a\ = - (Ba * Ca)3 (3 B\ + 3 B 2 C 2 - C\ B 2 - B\ C 2 ).Z640 (2.36) 

On the other hand, under the condition C\ = —3 Si, the condition (??) for the second 
possibility becomes 

C 2 Bf — 3C\BfB 2 — 3C 2 B\B\ + C\B\ = C 2 Bf + C\B\B 2 — 3C 2 B\B\ — 3B\B'2 

= - B 1 (3 Bl + 3 B 2 C 2 - C 2 B 2 - B\ C 2 ) = 0, 

which implies, by Eq. (??), 3 (a 4 + 2) (a 4 + l) 2 — (5 a 4 + 6) a| = for B\ ^ 0. Hence, 
according to Zoladek's classification (see (??)), this case should be included in the case. 
However, one can not prove this by directly using the conditions in (??) as well as that for 
the case (see Theorem 1.1). One must trace back to the original system coefficients. 

In [20], Zoladek used Bautin's system to verify his classification. Bautin's system is 
described by [6] 



dec 

— - = Ai x - y + A 3 x 2 + (2 A 2 + As) x y + A 6 y 2 , 

f mi 

— = x + Ai y + X 2 x 2 + (2 A 3 + A 4 ) x y - X 2 y 2 . 



(2.37) 



It is seen from (??) and (??) that Bautin's system has only 6 parameters, while Zoladek's 
system has 7 (in real domain) parameters. This indicates that Zoladek's system has one 
redundant parameter. In fact, putting Bautin's system in Zoladek's complex form gives the 
following expressions: 

A = Ai, A = U\ 3 + A 4 - A 6 -iA 5 ), S = -|(A 3 -A 6 ), 



C 



- (3A 3 + A 4 + A 6 ) +z(4 A 2 + A £ 



Then, applying the formulas given in (??) will immediately generate the centers conditions 
obtained by Bautin [6j. The above expressions clearly show that B 2 = 0. As a matter of 
factor, the integral factor for the system, corresponding to the second possibility, i.e., when 
A = A - 2 B = lm(B 3 C) = 0, is given by 



(B 2 x + B lV ) 



5S 1 (S^-3b|) 



C 1 (B'f+B^)-2B 1 (B'f-3B^) 



For B 2 = 0, the above expression is reduced to 

5 Bi 



1-2(2B 1 -C 1 )y 



1-2(A 1 -C\)y 



(due to Ay = 2 B x ) 



which is the integrating factor for the system, as shown in (??). 

Now we return to system (??). Among the four classifications, the Hamiltonian system 
(Qf) has been completely studied in [21], [22] : the system can have maximal two limit cycles. 
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In this paper, we will concentrate on the Q R - reversible case. Special cases for the reversible 
system have been investigated by a number of authors (e.g., see [23J, [2U [25j I2S1 [27]). It is 
easy to see that system (??) is invariant under the mapping (t, y) — » (—t, —y), where a\ 
and CI4 can be considered as perturbation parameters. The singular point (1,0) of (??) is 
a center when a\ < — 1; but a saddle point when a\ > — 1. a\ — — 1 gives a degenerate 
singular point at (1,0). Further, it is easy to verify that when {a\ + 1) > 0, there are no 
more singularity; while when {a\ + 1) 04 < 0, there exist additional two saddle points, given 
by 

_ J_ ± v /-a 4 (ai+l) \ 



(a ,2/ 



ai ai a4 



04 = is a critical value, yielding the two additional saddle points at infinity: (x* y*) = 
(— — , ±00). In summary, the distribution of singularity of the reversible system (??) has 
the following possibility (see Fig. [TJ where 1C+1S stands for one center and one saddle point, 
similar meaning applies to 2C, 2C+2S and 1C+3S): 

Two centers when a\ < — 1 and 04 < 0; 

Two centers and two saddle points when a± < — 1 and 04 > 0; 

Z643 (2.38) 

One center and one saddle point when a\> —1 and 04 > 0; 
One center and three saddle points when a\ > — 1 and 04 < 0. 

In this paper, we pay particular attention to ^ < - 1, a 4 < 0, for which system (??) 
has only two singularities at (0,0) and (1,0), both of them are centers. 

By adding quadratic perturbations to system (??) we obtain the following perturbed 
quadratic system: 

$= V (! + ai x) + e P(x, y) 

= y (1 + ax x) + e (a 10 x + a 01 y + a 2 o x 2 + a n xy + a 02 y 2 ) , 



= — x + x 2 + CI4 y 2 + e Q(x, y) 

= -x + x 2 + a 4 y 2 + e (b 10 x + b 01 y + b 20 x 2 + b n xy + b 02 y 2 ) 



WU (2.39) 



where < £ < 1, a^'s and fr^'s are perturbation parameters. 

Remark 1.4. The special system considered in [23] is the system (??) with 

a = — 3, c = — 2, b = 1. 

This is equivalent to our system when a\ = — 2 and = — 3 for which the system has 
only two centers at (0, 0) and (1, 0). Consider the 01-04 parameter plane, as shown in Fig. [TJ 
It can be seen that the case considered in [23] is just a point, (ai,a 4 ) = (—2, —3), in the 
parameter plane, marked by a blank circle in the third quadrant on the line 04 = | a% (see 
Fig. CD . 



P. Yu & M. Han 



4 limit cycles in near- quadratic systems 



15 



The special system studied in [23] is the system (??) with 

o = — 3, c = — 2, b — — 1. 

This is equivalent to our system when a\ = 2 and 04 = 3, for which the system has one 
center at (0,0) and one saddle point at (1,0). Thus, this case considered in [23] is again a 
point, (ai, 04) = (2, 3), in the 01-04 parameter plane, marked by another blank circle in the 
first quadrant on the line o 4 = | Oi (see Fig. [TJ). 

The cases considered in [251 [26] correspond to the system (??) with o = — 3, c = — 2, 
and be (-00, -1) U (-1, 0) in [25], and b G (0, 2) in [2B]. 

When e = in system (??), one can use the following transformation: 

x = v y = v 



to transform system (??) e=0 to 



dx / c 
*=»( 1+ 6 



- , ~2 . fl ~2 
rft = - X + X + b V -< 



Z645 (2.40) 



which is our system (??) with 
Equation (??) yields 



oi = 7, o 4 = T.IH6 (2.41) 


o 4 = -oi (6^0),/M7 (2.42) 

c 

which represents a line in the 01-04 parameter plane, passing through the origin with the slope 
-. In particular, the parameter values: a = — 3, c = — 2, 6g (—00, —1) U (—1,0) U (0,2), 
yielding oi = — | and 04 = — |, correspond to a part of the line, described by 

3 

O4 = — Oi Vai G (—00, —1) U (0, 00), lb48 (2.43) 

as shown in Fig. (TJ where the dotted line for Oi G [—1,0] is excluded from the studies [25], 26J . 

It should be noted that when a = — 3, c = — 2, the point (0, |) is a saddle point if and 
only if 

c 2 

l + - = l-->0 6 G (-oo,0) U (2, +00). 



Thus, the case considered in [25] has one center and one saddle point; while the case studied 
in [26] has two centers. But even these two studies together do not cover the whole line 
a 4 = I oi (the missing part is denoted by a dotted line segment in Fig. CQ . 

Another alternative form for a special case of our system (??) considered in [30] is 
described by 

dx r _ .„ . / l\n 

[l + 2 (1 - e) (x + -J 



x + dx 2 + ey 2 , 



lbA9 (2.44) 
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where e and d (7^ 0) are parameters. This system has a saddle point at the origin and a 
center at (x, y) = (— |, 0). Based on the two parameters, seven cases are classified |30j. We 
can apply the following transformation: 



* = i(*-i), V=\y, 



to system (??), yielding 



dx 



, 2(l-e) 
1 + — -x 

a 



d2/ , 2 , e 2 



' /650 (2.45) 



which has a center at the origin and a saddle point at (1,0). Then, setting 

2(~\ — p) e 
a x = K , ; , a 4 = -, Ml (2.46) 
a a 

in system (??) leads to our system (??). Equation (??) denotes a line, given by 

a 4 = r Oi,f&52 (2.47) 

2(1 -e) 

in the a\-a± parameter plane, passing through the origin with the slope 2 . However, it 
is easy to see that using our system (??) in analysis is simpler than using system (??). In 
fact, all the seven cases classified in [30] together denote a region in Fig. [fl see the shaded 
area in this figure. This area covers most of the region, defined by a\ > — 1. But the study 
given in [30J for the seven cases is restricted to local analysis on the bifurcation of limit cycles 
near a homoclinic loop, except the two lines (see Fig. [T]): 

a 4 = a 1 V01 G (-1,0) U (0,oo) ,1653 (2.48) 

which corresponds to the parameter value e = |, and 

a 4 = -^ai Vai e (0,oo),Z654 (2.49) 

which corresponds to e — > ±00. It has been shown [3D] that except the above two lines, for 
the parameter values in the shaded area, system (??) can have at most 2 limit cycles near a 
homoclinic loop under quadratic perturbation. 

Figure [H shows the a\-a± parameter plane associated with the reversible system (??), 
where the above mentioned case studies are indicated on the line a 4 = | a\ as well as in 
the shaded area. More precisely, a complete global analysis given in [25J, which includes the 
result in [23] as a special case, shows that corresponding to each point on the line segment 
a 4 = § a,\ (ai > 0), the system has one center and one saddle point, and has maximal 2 limit 
cycles. In [30] it is shown for each point in the shaded area (except the two line segments 
a 4 = a\ (ai > —1) and a 4 = —\cl\ (a\ > 0)), which contains the above line segment, 
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a 4 = |(ai — 5) 




0,4 — — ~ CLi 



0.4 = |(6ai + 5) 

Figure 1: Case studies for the Qf reversible system. 



the system has one center and one (or three) saddle(s), and has maximal 2 limit cycles, 
but restricted to local analysis near one homoclinic loop. Similarly, a global analysis given 
in [26], which contains the result in [23] as a special case, proves that corresponding to each 
point on the line segment — |a x (oi < — 1), the system has two centers, and exhibits 
maximal 3 limit cycles around one center. The technique of Poincare transformation and 
Picar-Puchs equation, used for the above mentioned global analysis on parameter unfolding, 
seems not possible to be generalized to consider general situation for arbitrary points in 
the 0,1-0,4 parameter plane. The two particular dash-dotted lines: 04 = | (a± — 5) Vc&j G 
(— 00, —1) U (—1, 00), and 04 = I (6ai + 5) W a\ G (—00, —1), as well as the five dark circles 
correspond to our results, presented in the next two sections. In particular, we will show 
that there exist 3 small limit cycles on the two dash-dotted lines, and at least 4 limit cycles 
for the parameter values marked by the five dark circles. 

In the following, we will use the perturbed quadratic system (??) for our study on 
bifurcation of limit cycles. Without loss of generality, we may assume (e.g., see [23]) that 
aoi = &20 = a n = a 02 = bio = &20 = bo2 = 0. Thus, system (??) is reduced to 



dx 

— = y (1 + ax x) + e aw x, 



%= -x + x 2 + a 4 y 2 + e(b iy + b n xy), 



/655 



(2.50) 



where a% < — 1 and 0<eCl. 
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3 Hopf bifurcation associated with the two centers 

In this section, we study Hopf bifurcation of system (??) from two centers (0,0) and (1,0), 
leading to bifurcation of multiple limit cycles. The result is summarized in the following 
theorem. 

Theorem 2.1. When a\ < — 1, the quadratic near-integrable system (??) can have small 
limit cycles bifurcating from the two centers (0,0) and (1,0) with distributions: (3,0), (0,3), 
(2,0), (0,2) and (1,1). (2,1)- or (1,2) -distribution does not exist. 

Proof. Consider system (??) for ai < — 1. The system (??) e=0 is a reversible integrable 
system. In order to compute the Melnikov function near the two centers (0, 0) and (1, 0), we 
need transform system (??) e= o to a Hamiltonian system. The integrating factor 7 is given 
in (??). Now, introducing dt = 7 dr into (??) yields the perturbed Hamiltonian system: 

% , , im (3 - 1} 

_ = j(-x + x 2 + a 4 y 2 ) + £7(f> i V + hi xy), 
dr 

with the Hamiltonian of (??) £ =o, given by 

T + ax — 04) (1 + 2 a 4 x) x 2 



H(x, y) = - sign(l + a\x) |1 + a\x\ 



y 2 + 



Z657 



04 (ai — a 4 ) (ai — 2 04) a\ — 

(3.2) 

for a 4 7^ 0, ai 7^ 04, a\ 7^ 2 04. The cases 04 = 0, a± = 04 or a± = 2a 4 will not be 
considered in this paper. 



Note that 

h 00 = H(0,0)= l + ai-a 4 for 1 + aiX>0; 

2 a 4 [ai — 04) (ai — 104) 

h 10 = H(i,o) = - ^ + 1 )( Q ; + 1 ) ( _i_ Ql) -^ i for i + aiX < . 

2 04 (ai — a 4 ) [ai — 2a 4 ) 



Z658 (3.3) 



Since in this paper, we concentrate on the case that system (??) e=0 has only two centers, we 
assume a± < — 1, 04 < 0. Thus, 

lim H(x, y) = + 00 and lim H(x, y) = — 00. 

X^- — ~ — + 

It is easy to see from system (??) that the trajectories of (??) £ =o rotate around the center 
(0,0) in the clock-wise direction, while rotate around the center (1,0) in the counter clock- 
wise direction, as shown in Fig. Thus, the values of h in H(x, y) — h are taken from the 
two intervals: h G (h 00 , 00) for 1 + a x x > 0, and h G (—00, h w ) for 1 + ai x < 0. It should 
be noted that is not necessarily larger than h w . The analyses on the two half-plane in 
the x-y plane (see Fig. [2]), divided by the singular line 1 + a± x — 0, are independent. 
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x 

Figure 2: A phase portrait of the reversible system (??) with two centers for a\ = —3, 

_ 8 
3' 



Next, introduce 

{h £ (hm, oo), for 1 + ai x > 0, 
h G (— oo, hxo), for 1 + a\ x < 0, 

and define the Melnikov function: 

M(/i, dij, kj) = f q(x, y, by) dx - p(x, y, cfy, Z&60 



(3.4) 



(3.5) 



where p(x,y,a,ij) = r yaiox and q(x,y,bij) = 7 (&01+&11 x) y. Using the results in [T9 | l3T | 132] . 
we can expand M near h = h 00 and h = h 10 as 



M (h, dij, bij) = fi 00 (h - h 00 ) + fa! (h - h 00 ) 2 + ^ 02 (h - h 00 ) 3 

+ /i 03 (h - /i 00 ) 4 + 0((h - /i 00 ) 5 ), for < h - h 00 « 1, 

Mi(/i, ay, %) = /iio (/iio - /i) + (/iio - /i) 2 + M12 (ho - /i) 3 

+ A*i3 {h 10 - hY + O((/ii - h) 5 ), for < h 10 - h < 1, 



mi (3.6) 



where the coefficients /x^, z = 0, 1; j = 0, 1, 2, • • • can be obtained by using the Maple 
programs developed in [33] as follows: 



^00 = 2 7r(ai + 6 i), 
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(1540 - 980ai - 280a 4 + 861a? - 1512aia 4 - 3948^ - 626a? + 1566a?a 4 
+1620aia^-4432a^+313a^-1018a?a 4 -279a?a^ + 3080aia^-2096a^)a 10 
+ (1540 + 700a! + 3080a 4 + 21a? + 168aia 4 + 2772a 2 - 2a? + 126a?a 4 
- 828a ia 2 + 1424a 3 + a\ - 58a? a 4 + 369a? a 2 - !Yla x a\ + 400a^) b 01 
+ 24 b u (1 + a 4 ) (70 + 35ai + 70a 4 + a? - 17 + 52a 4 ) 6 n ] , 

(3403400-300300a 1 + 3003000a 4 +690690a?-4984980aia 4 -7327320a 

- 500885a? + 3314850a?a 4 - 4430580aia^ - 17811640a^ + 323121a^ 
-2444439a?a 4 +4201218a?a 2 + 5794692a 1 a 4 i -18033936a 4 ! -168603a? 
+ 1420500a?a 4 - 3253551a?a 4 : - 1296282a?a 3 + 121079040^ 

- 10462368a 4 ' + 56201a? - 52031 la?a 4 + 14712870^ - 407053a? a 3 

- 4589772a?a^ + 7149264a 1 a 4 ' - 3159616a^) a w 

+ (3403400 + 3303300a! + 10210200a 4 + 690690a? + 5825820aia 4 
+14294280a 4 ; + 11935a?+404250a?a 4 + 2721180a 1 a 4 ; + 12236840a 4 i 

-699a^ - 11379a?a 4 + 262458a?a 4 : - 1891308a!a^ + 6994704a^ 

+417a? + 1380a^a 4 - 149091a?a^ + 1121838a?a 4 : - 29645760^ 

+2670432a 4 ' - 139a? - 291a?a 4 + 46227a*a 2 - 366193a?a 4 i 

+1076988a?a 4 t - 1335216a 1 a 4 ' + 578624a! 1 ) b 01 

+ (3603600 + 3603600ai + 10810800a 4 + 790020a? + 6597360aia 4 

+15024240a 2 + 12600a?+480060a?a 4 + 3764880a 1 a 2 + 12514320a 4 i 



864 



7T 



622080 



+180ai - 10800a?a 4 + 11340a?a4 - 6184800^ + 6566400a 



4 



+180a^a 4 -23400a?a^ + 321300a?a^-1389600aia^ + 1869120a2) b u 



2 TT (-1 - ai ) 3/2 [(1 - 2 04) aio + (1 + ai) (6 i + 611), 



7T _ 2 ( a i~"4) 



x 



(10 + 33ai — 6a 4 + 36a? — 21aia 4 — 24a?a 4 + 30aia 2 — 804) a w 
+ (1 + ai) (10 + 21ai - 10a 4 + 12a? - 15aia 4 + 4a 2 ) h 



01 



-(l + ax)(l 



(2 + 3ai - 4a 4 ) b n 



7J- _ (5a 1 -8a 4 ) 

864 ( - 1 -°' r 



X 



(1540 + 7140ax - 2800a 4 + 13041a? - 11592aia 4 + 2212a 2 + 11448a? 
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-18072^04 + 86280104 - 1 112o|+ 752a 4 - 12024a?a 4 + 12213a?a| - 5232aia^ 
+4320a? - 1728a?a 4 + 6192a?a^ - 7938a?a^ + 4272aia* - 800a^) a w 
+ (1+ai) (1540 + 5460ai-3080a 4 + 7161a?-9072a 1 a 4 + +2772a2+4104a? 

- 9030a?a 4 + 6372aia^ - 1424a^ + 864a? - 3096a?a 4 + 3969a?a;[ 

- 2136aia^ + 400a^) b m 

- (l + ai)(l + a 4 ) (140 + 420ai - 420a 4 + 423a? - 996aia 4 + 576a^ 
+ 144a? - 633a?a 4 + 8880^ - 400a|) b u 



7f _ (2 (a 1 -3a 4 ) 

^ u 1244160 V ' 

x "(3403400 + 20720700ai-9809800a 4 + 53243190a?-54234180aia 4 + 13093080a^ 
+ 74334645a? - 123735150a?a 4 + 65571660aia^ - 10776920a^ + 60023916a? 

-147900519a?a 4 + 131934978a?a^-49682268aia^ + 6439744a^ + 27002160a? 
-95460120a?a 4 + 132380865a?a^-89408610a?a^ + 29027880aia^-3527040af 
+ 5443200a? - 28946160a?a 4 + 63998532a?a^ - 74879613a?a^ +48498336a?a^ 
- 162963360^ + 2181248a^ - 1555200a?a 4 + 9603360a?a^ - 24061752a?a^ 
+ 31232358a?a^ - 22072536a?a^ + 8011296aia^ - 1157248a^) a w 
+ (1 + ai) (3403400 + 17117100a 1 -10210200a 4 + 35225190a?-45225180aia 4 
+14294280a^ + 37785825a? - 79202970a?a 4 + 54455940aia^ 
-12236840a^ + 22125636a? - 68371209a?a 4 + 77864598a?a^ 
-38601828aia^ + 6994704a^ + 6629040a? - 28984608a?a 4 
+49687587a?a^ - 41614974a?a^ + 16953984aia| - 2670432a^ 
+777600a? - 4801680a?a 4 + 12030876a?a^ - 15616179a?a^ 
+11036268a?a^ - 4005648aia| + 578624a^) b 01 
- (l+ai)(l+a 4 ) (200200 + 900900ai-800800a 4 + 1600830a?-3132360aia 4 
+ 15307600^ + 1397655a? - 4596480a?a 4 + 50085000^^ 
-1808240a^ + 594864a? - 3001266a?a 4 + 5594022a?a^ 
-4568112aia^ + 1379936a^ + 97200a? - 736776a?a 4 
+2162079a?a2-3080268a?a^ + 2136528aia^-578624a^) b u 



Remark 2.2. The coefficients fi^j listed above are applicable as long as (0, 0) is a center, 
and the coefficients are applicable as long as (1,0) is a center, regardless the number 
and distribution of the system's singularities. Therefore, for each point on the whole line 
a 4 = ^ (ai — 5) (see Fig. [I]), there always exist 3 small limit cycles bifurcating from the center 
(0,0), no matter whether the system has two centers, or one center and three saddle points, 
or one center and one saddle point. For each point on the line segment a 4 = | (6 ai + 5) {a\ < 
— 1), the system can have 3 limit cycles bifurcating from the center (1,0). This indicates 
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that the results given in [23l |24~ , 30j showing that the reversible near-integrable systems with 
one center and one saddle point can have maximal 2 limit cycles is conservative, since on 
the part of the line = ~ (ai — 5) in the first quadrant (ai > 5) such a system can have at 
least 3 limit cycles. 

First, we consider the maximal number of limit cycles which can bifurcate from the 
center (0,0). Setting yU o = yields 

b i = — a 10 ,lb62 (3.7) 

and then we have 

/An =71" («i — 1 — 04) (01 + 204)010 + (1 + (14) 6n J663 (3.8) 
In order to have /i i = 0, we suppose a 4 7^ — 1 and choose 

(ai — 1 — a A ) (ai + 2a 4 ) 



ai ./664 (3.9) 



/665 



on :— 

1 + 04 

Then, /i 2 and /i 3 are simplified to 

//02 = — ai (ai — 04) (ai + 2 04) (ai — 3 04 — 5) aio, 
3 

= ~ Tii ai _ a4 ^ + 2a ^ ^ 770 + 105ai + 1400a4 + 42a i ~ 434a l fl 4 

+1274al - 13a? + 128a?a 4 - 415aia^ + 444a^) ai . 

(3.10) 

There are five choices for ^2 = 0. Except the choice CL\ — 3 (24 — 5 = 0, all other choices lead 
to fi 0i = 0, i = 3, 4, • • •. Thus, letting 

a 4 = - (ai - 5),/666 (3.11) 
3 

which implies a x 7^ 2 when 04 7^ — 1. Since we assume ax < — 1, for this case (i.e., when 
the condition (??) holds), 7^ — 1 is guaranteed. Then, we have 

25tt 

^03 = -T77T fl i ( a i + !) ( a i ~ 2 ) 2 ( 2 a i + 5 ) a io, 
lo2 

^ = K + l)(ai-2) 2 (2 ai + 5) (oi+4)(17ai + 518)oi 

o 1 4o 

£*io = - (-1 - «i) _3/2 ai (2 ai + 5) a w , 

257T _ 2 ( 2 °i+ 5 ) 

A*n = ( _1 _ ai ) 3ai a i ( a i - 2 ) 2 (2 «i + 5) aio, 



implying that in addition we need 
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(2ai + 5)ai + 0./667 (3.12) 

Under the above conditions (??), (??), (??) and (??), we obtain /x o = A*oi = /^02 = 0, 
but /i 03 7^ 0, /xio 7^ 0. Hence, at most 3 small limit cycles can bifurcate from the center (0.0) 
with no limit cycles bifurcating from the center (1,0). Further, giving proper perturbations 
to the parameters a 4 (or ai), bu and &oi, we can obtain 3 small limit cycles bifurcating from 
the origin. This shows that the conclusion is true for the case of (3, 0)-distribution. 

Next, consider the (0, 3)-distribution. Similarly, letting fi w = yields 

boi = - fen + ~ 1 a 10 ./668 (3.13) 

1 + Oi 

Then, /i n becomes 

A*n = 7r(-l-ai)~ ( "i 4) (ai + 2a 4 )(2ai-a 4 + l) ai -(l + ai) 2 (ai-a4 + l) &n -^69 (3.14) 
Hence, we set 

7 (ai + 2a 4 )(2ai - a 4 + 1) , . , , 

t>n = — — 777 — — aio, (ai - a 4 + 1 7^ 0),/670 (3.15) 

(1 + ai) 1 (ai - a 4 + 1) 

to yield // n = 0, and 

A*i2 = - (-1 - 01) 2ai ai (ai - a 4 ) (a 4 + 2 a 4 )(6 01 - 3 a 4 + 5) a 10 , 

/ii 3 = — (-1 - ai)~ 2( V a4) ai (a 1 -a 4 )(ai + 2a 4 ) (770 + 2205ai - 1400a 4 + 2142a? lb71 
288 

- 3234aia 4 + 1274af - 720a? - 1962a 2 a 4 + 1689oial - AAAa\) a 10 . 

(3.16) 

The only choice for n\ 2 = is 6ai — 3a 4 + 5 = 0, from which we have 

a 4 = ^(6ai + 5).Z672 (3.17) 
3 

This implies that ai — a 4 + 1 = — (ai + f ) > for a\ < — 1. Further, we obtain 

257T 10+11 ai 

/X13 = - — (-1 - ai) "i 01 (3ai + 2) 2 (3ai + 5) aio, 

5"7T 80+87ai 

A*i4 = -tt^ (- 1 - fl i) 601 ai(3ai + 2) 2 (3ai + 5)(3ai + 4)(501ai + 518)ai 
1749o 

10tt 

^00 — o /1 1 \9 fll 0l + ^) ai 0' 

3 (1 + ai)" 1 
257]- 

V01 = ~ 324(1 + ai )2 fll ( 3 ai + 5 ) ( 3 ai + 2 ) 2 ai0 ' 
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implying that in addition we require 

(3ai + 5)ai ^ 0./o73 (3.18) 

Under the above conditions (??), (??), (??) and (??), we have /iio = /in = M12 = 0, 
but yU 13 7^ 0, Moo 7^ 0. Further, by properly perturbing the parameters a 4 (or ai), b n and 
601, we can obtain 3 small limit cycles bifurcating from the center (1,0), but no limit cycles 
from the origin. This proves the case of (0, 3)-distribution. 

For the case of (2, 0)-distribution, it follows from the conditions (??) and (??), and 
a 4 7^ — 1 that Moo = Moi — 0, and 

7T 

M02 — — cli (a>i — a 4 ) (ai + 2 a 4 )(ai — 3 04 — 5) aio, 
3 

2 7T 

M10 = - 7— — 77—; vTT? fll ( Gl ~~ ° 4 ) ( Gl + 2 ° 4 ) ai °- 

(1 + a 4 J (-1 - aij 13 / 2 

Thus, M02 7^ implies M10 7^ 0, indicating that the conclusion holds for the case of (2,0)- 
distribution. if 04 = — 1. 

When a 4 = — 1, (??) becomes 

Moi = 7r ai (ai — 2) aio 7^ for ai < — 1 and a±o 7^ 0. 

Under the conditions o i — — a io an d a 4 = — 1, M10 and fj, n becomes 



A*io = - 2 ?r (-1 - 01) 3/2 [(ai - 2) aio - (1 + ai) 611 
M11 = 7T (-1 - ai) "i ai(ai-2)ai , 



/674 (3.19) 



which shows that Mil 7^ for ai < — 1 and a 10 7^ 0. But we can choose 

ai-2 

011 = 7-7 aio 

1 + ai 

to obtain /x 10 = 0. Thus, for this case we have a (1, l)-distribution. 

Similarly, for the (0, 2)-distribution, we use the conditions (??) and (??) to obtain 

M12 = - (-1 - «i) 2ai 01 (ai - a 4 ) (ai + 2a 4 )(6a! - 3 a 4 + 5) ai , 

2 7T 

Moo = 77- 777 77-7 ai (ai - a 4 ) (a : + 2 a 4 ) ai . 

(1 + ai) z (ai - a 4 + 1) 

This indicates that M12 7^ implies Moo 7^ 0, and so the conclusion for the case of (0, 2)- 
distribution is also true if a\ — a 4 + 1 7^ 0. 

When ai — a 4 + 1 = 0, i.e., a 4 = a x + 1 < 0, (??) is reduced to 

2 (a^ — a-4) 

M11 = ^l - 1 ~~ a i) a i(3 o-i + 2) a 10 7^ for ai < — 1 and ai 7^ 0, 
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and /too and /toi become 

2 TT 

-[(ai + 2)aio-(l + ai)6nJ, 

W75 (3.20) 



A*oo 



1 + ai 



7T 

A*oi = - z «i (3 ai + 2) aio, 

1 + ai 

which clearly shows that /x i 7^ for ai < — 1 and a 10 7^ 0. However, we may choose 

a x + 2 

on = — aio 

1 + ai 

to obtain // o — 0. Thus, for ai — a 4 + 1 = 0, we have a (1, 1) -distribution. 

Finally, suppose the condition given in (??) is satisfied, i.e., 6 01 = — aio, then substitut- 
ing this into fiw to solve bu to obtain 

o n = ^±^i./676 (3.21) 
1 + ai 

Then, under the conditions (??) and (??), we obtain 

TT 

(Mn = z— t ai (ai - a 4 ) (a : + 2 a 4 ) ai , 

1 + a i , /&77 (3.22) 

_ 2(a 1 -a 4 ) V / 

^11 = - 7r (-1 - ai) ai ai (ai - a 4 ) (ai + 2a 4 ) aio, 

which shows that /x i 7^ implies /in 7^ 0, and thus in general the conclusion is true for the 
case of (1, l)-distribution. 

As we have seen in the above analysis, if the condition (??), a 4 = | (a x — 5), is not 
used, then we can only have 2 limit cycles bifurcating from the origin, but no limit cycles 
can occur from the center (1,0). In other words, we can obtain one more limit cycle, by 
using the condition a 4 = | (a x — 5), only bifurcating from the center (0,0). Similarly, if the 
condition (??), a 4 = | (6a x + 5), is not used, then we can have only 2 limit cycles bifurcating 
from the center (1,0), but no limit cycles can bifurcate from the origin. Then, condition 
a 4 = I (6ai + 5) can be only used to get one more limit cycle around the center (1,0), rather 
than the origin. Therefore, (2, 1)- or (1, 2)-distribution is not possible. 

This completes the proof of Theorem 2.1. □ 



4 Limit cycles bifurcating from closed orbits 

In this section, based on the results of the small limit cycles obtained in the previous section, 
we want to investigate the possibility of existence of large limit cycles by applying the 
Melnikov function, defined in (??). We have the following result. 

Theorem 4.1. For the case of bifurcation of small limit cycles from the two centers (0, 0) 
and (1,0) with (3,0) -distribution (respectively, (0,3) -distribution) there exists at least one 



P. Yu & M. Han 



4 limit cycles in near- quadratic systems 



26 



large limit cycle near Lh for some h G (—00, hio) (respectively for some h G (hoo, 00) ). For 
the case of limit cycles with (2 , 0) - distribution (respectively, (0,2) -distribution) there exist at 
least two large limit cycles, one near L hl for some hi G (—00, h w ) and one near L h2 for some 
h-2 G (/loo, 00). The corresponding values of the parameters a± and a 4 for the existence of A 
limit cycles can appear at least in some regions in the 01-04 parameter plane. 

Remark 4.2. Theorem 4.1 gives a positive answer to the open question of existence of 
limit cycles in near-integrable quadratic systems: at least 4 limit cycles can exist. For the 
case of (1, 1) -distribution, so far no more large limit cycles have been found. 

Proof. It follows from (??) with 

0-^+20.4 aj+2a4 

p(x,y,a,ij) = |1 + aix\ ~ a w x, q(x,y,b ij ) = |1 + ai x\ ~ (b 01 + b n x) y 



that 



M(h, ai, 04,010,601,611) 

f q(x,y,bij)dx - p(x,y,aij)dy 

J L h 

f q(x,y,b i:j )dx - f p(x,y,a i:j )dy 

JL h JL h 

f q(x,y,b i:j )dx+ f yp x (x,y,a i:j )dx 



q(x,y,b i:j ) +yp x (x,y,a i:j ) 



dx 



01+204 01+204 / (a 1 + 2a 4 )x" 

|l + aix| a i (6 i + 6n x) + |1 + ai x\ a i aio 1 =F 



ai+2a.4 

|1 + ai x\ 



|1 + Oi x\ 

(a w + 6 i) + 611 x — sign(l + ai x) a w (01 + 2 a 4 ) 



y dx 



x 



1 + Oi x\ 



y dx 



= (aio + 601) Io(h, ai, a 4 ) + 6n h(h, a 1 , a 4 ) + ai h{h, ai, a 4 ) 

M (h, ai, a 4 , aio, b i, 6n) for og(oqo,oo), when 1 + ai x > 0, 



Mi{h, ai, a 4 , aio, 6oi) 611) for /i G (—00, /iio), when 1 + a! x < 0, 



/678 



(4.1) 



where 



= < 



r 01+204 
/ (/i, ai, a 4 ) = a> |l + Oix| "i yda; 

lax _ 01+204 

(1 + ai x) °i y+dx, V A G (/ioO) 00)) when 1 + OiX > 0, 

n 

—2 / (— 1 — aix) 2/ + dx, V /i G (— oc, /^io), when 1 + < 0; 



h(h, ai,a 4 ) = f 



|1 + ai x| ai xy dx 
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[1 + aix) °i xy+dx, V h G (/too, oo), when 1+aix > 0, 
^— 1 — aix) a i xy+dx. V /i G (— oo, /iio), when 1 + aix < 0; 



/■ 2(a 1 + a 4 ) 

J 2 (^, ai, a 4 ) = — (ai + 2 a 4 ) sign(l + aix) ® |l + aix| a i xydx 



-2(ai + 2a4 



2( ai + a 4 ) 



^1 + aix) a i xy+dx, V/iG (/iqO) oo), when l + aix>0, 



2(a 1 +a4) 



2(a 1 + 2a 4 )/ (— 1 — aix) a i xy+dx, V/tG (— oo, /iio), when l + aix<0. 



Here, 



2/+ 



(1 + ai — 04) (1 + 2 a 4 x) 
a 4 — a 4 a 4 (ai — 04) (ai — 2 a 4 ) 



4t 

+ 2h sign(l + ai x)\l + ai x\ a i 



1/2 



W79 (4.2) 



and x m i n and x max are solved from the equation, y+ = 0, for /i G (/ioo, 00) when 1+ai x > 0, 
and for /i G (—00, /iio) when 1 + a 4 x < 0. 

Since one can not find the closed form of the integrals Ii(h, 04), i = 0, 1, 2, for 
general a\ and a 4 , nor the technique of Picard-Fuchs equation can be applied here, we shall 
choose some values for a\ and a 4 and then find numerical values of the integral. We first 
use the results given in the previous section to determine b i, b u , and a 4 , and then choose 
proper values for a\ to find more limit cycles. 

(A) First, consider the (3, 0)-distribution. For this case, we have 



a4 = ~(ai-5), 6 i 



aio, &11 = - 10 (1 + 01) a 10 . 

which denotes a point (a blank circle) on the line 
a 4 = I (ai — 5) in the a 4 -a 4 parameter plane (see Fig. [p. Further, we have 6 n = y a 10 , 
and 

30 N-f , ,7 



Taking a± = — y yields a 4 



65 
21' 



7=(l-y 
Then, the Hamiltonian (??) becomes 



x 



(x + —). 
V T 30 ; 



rr/ , 16807 (16250 y 2 + 13650 x 2 + 2730 x- 441) r . 7 

xii ( x t/ ) = tor X =P — 

K ,yj 32500 (7- 30 x) (40353607- 172944030 x) 4 / 9 r 30 ! 



with 



441 



00 



> h 



33957 /23x5/9 



32500 ' '" 10 747500^ 7 
The Melnikov functions Mi(h,aio) can be expressed as 



Mi(h, a w ) = M i0 (h) a 10 , i = 0, 1.Z680 



(4.3) 
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Figure 3: Functions M 00 (h) and M w {h) under the conditions /ioo = /An = H02 - 
and ^10 7^ 0, for a± — — y and 04 = |(ai — 5) = — 1|: (a) M 00 (h) > for ft, G [/too 
^0 = - 3U0 ~ "0.01357; and (b) M 10 (ft) for ft G (-00, h], with ft 10 = - |§| v . 
-0.08797, crossing the ft-axis at ft = ft? G (-0.9250363254, -0.9250363253). 



0, /i03 7^ 
+00), with 

23N5/9 „ 
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Figure 4: Illustration of the existence of 4 limit cycles when a± — — y, = |(ai — 5) = 
-|f - £i, and do = 2^0, b u = ™ a 10 - e 2 , b 01 = - a 10 - e 3 , where < e 3 < e 2 < e x < e: 
(a) An unstable large limit cycle enclosing the center (1,0); and (b) Zoomed area around 
the center (0,0) showing the existence of 3 small limit cycles. 
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Without loss of generality, we may assume 

a w >0,lb81 (4.4) 

and thus Mi(h,a w ) and M i0 (h) have the same sign. It is noted that for the above chosen 
parameter values, we have 



139150000 7T n , 2500 v / T6Tvr 

am = am > and uin = Qin < 0. 

P 453789 P 3703 

The computation results of M o(h) for h G (/ioo,oo) and Mi (/i) for h G (— oo, hio) 
are shown, respectively, in Figs. |3]^a) and [3](b) . Figure [3]^a) shows that M 00 (/i) > for 
h G (/loch oo), and its sign agrees with that of fio3 > for < h — h 00 <C 1, as expected. It 
is also noted, as shown in Fig. El^b), that the sign of M w (h) agrees with that of /jLiq < for 
< h w — h <C 1. However, unlike the interval h G {h 0Q , oo), this interval contains a critical 
value h = h\e (-0.9250363254, -0.9250363253) at which M w (h\) = and the function 
M 10 (h) changes its sign as h crosses this critical point. Thus, for this case, besides the 3 
small limit cycles, there exists at least one large limit cycle bifurcating from the closed orbit 
Lh* of (??)• This large limit cycle is shown in Fig. ID^a), which encloses the center (1, 0); and 
Fig. IHb) illustrates the existence of 3 small limit cycles around the center (0, 0). 

(B) For the case of the (0, 3)-distribution, we have 



1 2a 4 -l (ai + 2a 4 )(2ai - a 4 + 1) 

-(6ai + 5), b 01 = -b 11 + — oio, b n = —— — — — 

3 1 + a\ (1 + ai) z [a\ — a 4 + 1) 



By choosing a\ = — ||, we have a 4 = — ||, 6 i = — a w and bn = a w- The point 
(ai, a 4 ) = (— ||, — || ) is marked by a blank circle on the line a 4 = | (6 a x + 5) in the ai-a 4 
parameter plane (see Fig. [1]). Moreover, 

7=(l-^)'- (s^), 

and the Hamiltonian is 

TT . . 345025251 (2750 y 2 + 9350 x 2 - 16830 x + 7803) r 51 

Jj ( x Vi — for x =?= — 

v ' y; 5500 (51 - 70 x) (897410677851 - 1231740146070 x) A ' 7 r 70 ' 

with 

7803 , 44217 / 19x3/7 

hnn = > Alio = — 

5500 104500^51 ' 

For this case, /ioo and /ii3 become 

10500 7T n , 4561235000 /51s 2/7 

Moo = Oio < and Ui 3 = — 71 a w > 0. 

p 361 p 565036352721 U9/ 

The computation results of M 00 (h) for /i G (/i 0jOo) and M w (h) for /i G (— 00, h\o) 
are shown in Figs. |5la) and|5]^b), respectively. As shown in Fig. [Sfa), the sign of Mqq(K) 
agrees with that of ^too < for < h — h 00 <C 1, and in addition the function M 00 (h) 
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150 




-5 -4 -3 -2 -1 

h 

Figure 5: Functions Mqq{K) and M w (h) under the conditions /iio = fin = Hn = 0, ^13 7^ 
and yUoo 7^ 0, for a\ — — |^ and 04 = |(6ai + 5) = — 1| : (a) Mqq(K) for /i G [/io, +00), with 
fr 00 = » 1.41873, crossing the fr-axis at h = h* 2 e (13.3847179116, 13.3847179117); and 
(b) M 10 (h) > for h G (-00, h], with h w = - ^Ci) 3/7 ~ - 0.27714. 
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Figure 6: Illustration of the existence of 4 limit cycles when a\ — — «4 = §(6ai + 5) = 
-§f-£i, andaio = 10, bu = ^ai -e 2 , b 01 = -^-aio + £ 3 , where < e 3 < e 2 < S\ < e: 
(a) An unstable large limit cycle enclosing the center (0,0); and (b) Zoomed area around 
the center (1,0) showing the existence of 3 small limit cycles. 
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crosses a critical value at h = h\ G (13.3847179116, 13.3847179117), at which it changes 
sign. Figure [5(b) shows that M w (h) > for h G (— oo,/tio), and its sign agrees with that 
of /ii3 > for < h w -/i<l. Hence, for this case, in addition to the 3 small limit cycles, 
there also exists at least one large limit cycle bifurcating from the closed orbit of (??). 
This large limit cycle is depicted in Fig. Efa), which encloses the center (0, 0); and Fig. E(b) 
illustrates the existence of 3 small limit cycles around the center (1,0). 

(C) Now consider the (2, 0)-distribution. For this case, the condition a 4 = | (ai — 5) is 
not used. We need to determine the values for both a\ and a 4 . We choose 

«i = -4, a 4 = -— , 
5 

which represents a point in the third quadrant of the a\-a± parameter plane (see the dark 
circle in Fig. [T]near the line a 4 = |(ai — 5)). Thus, 

14 



7 =(l-4*)~ 1;f { X ±h 



In addition, we have b m = — a w , bu = ^ a w , and 

„, x 192 y 2 + 480 x 2 - 180 x + 25 r ,1 

H(x,v) — : tt-tz for x 4" ~, 

v ,yj 384(1 -4x) 9 / 5 T 4' 

with 

t 25 325 i /r 

h 00 = > h 10 = 3 1/5 . 

00 384 3456 

For this case, /i 02 and fi w are reduced to 

1344 40^ 
^ m = " ~125 71 aw < = ^— 7raio<0. 

The computation results of M 00 (h) for /i G (/ioo,oo) and Mi (/i) for h G (— oo,/iio) are 
shown, respectively, in Figs. [7(a) and[7Jb). As shown in Fig. UJa), the sign of M 00 (h) agrees 
with that of /X02 < for < h — h 00 1. Moreover, the function M 0Q (h) crosses a critical 
value at h = h% G (0.1448192224, 0.1448192225) at which it changes sign. Figure[7jb) shows 
M w (h) for h G (— oo,h w ), whose sign agrees with that of fi\o < for < h w -/i< 1. 
Also, M 10 (h) crosses a critical value at h = hi G (-0.5822537644, -0.5822537643) at 
which it changes sign. Therefore, for this case, besides the two small limit cycles, there exist 
at least two large limit cycles bifurcating from the two different closed orbits Lh* and of 
(??). One large limit cycle surrounding the center (1, 0) is shown in Fig. EJa), while another 
large limit cycle enclosing the center (0, 0) with 2 small limit cycles is depicted in Fig. E^b). 

(D) Finally, consider the (0, 2)-distribution. For this case, the condition a 4 = | (6 ai+5) 
is not used. Taking 

4 6 

^1 — — - , 04 — — — , 
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-1.5 I ' ' ' 1 

-0.8 -0.6 -0.4 -0.2 

h 

Figure 7: Functions M 00 (h) and M w (h) under the conditions //oo — A*oi ^02 7^ and 
fi W 7^ 0, for cii = — 4 and a 4 = — y: (a) M 00 (h) for h G [/ioo, +00), with ^0 = H; ~ 
0.06510, crossing the /i-axis at h = h* 3 G (0.1448192224, 0.1448192225); and (b) Mi (/i) 
for /i G (—00, /ii], with /i 10 = — Hf^ 1 / 5 « —0.11715, crossing the /i-axis at /i = /14 G 
(-0.5822537644,-0.5822537643). 
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Figure 8: Illustration of the existence of 4 limit cycles when a\ — — 4, 04 = — ^, and 
aio = ^11 = W ai ° ~ £l ' anc ^ ^ 01 = — a io — £ 2, where < 62 £i s: (a) An unstable 
large limit cycle enclosing the center (1,0); and (b) Zoomed area around the center (0,0) 
showing the existence of 1 large limit cycle and 2 small limit cycles. 
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Figure 9: Functions M 00 (h) and M"i (/i) under the conditions fj,io = /in =, ^ and 
yUoo 7^ 0, for ai = — | and 04 = — |: (a) M 00 (h) for fa G [/too, +00), with /i = f§§ ~ 
2.53096, crossing the /i-axis at /i = hi G (12.6197809949, 12.6197809950); and (b) M 10 (/V) 
for /i G (—00, hi], with /i 10 = — 3 4 / 5 ~ —1.41107, crossing the /i-axis at /i 
(-3.1388150376, -3.1388150375). 
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Figure 10: Illustration of the existence of 4 limit cycles when a± — — g, a& = — |, and 
a w = 1, &n = aio — £i, &oi = _ ^§ a w + £2, where < £2 £1 «C £: (a) An unstable 
large limit cycle enclosing the center (0,0); and (b) Zoomed area around the center (1,0) 
showing the existence of 1 large limit cycle and 2 small limit cycles. 
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yields 

7 = ( 1_ 3 X ) 5 ( x ^l)- 

The point (— |, — |) is marked by a dark circle near the line a± = | (6ai + 5) in the ai-a4 
parameter plane (see Fig. [p. Further, we have o i — — a io 5 ^11 — ^nr a icb an d 



rr/ . 243 (64 y 2 + 480 x 2 - 780 x + 325 r ,3 
H(x,y) = — — — — tor 

v ' y; ((3 -4 a;) (324 a; - 243) 4 / 5 r 4 

with 

° u 128 128 
For this case, fi 00 and /ii 2 are simplified as 

896 „ , 448 



^00 



7r aio < and fi u = - — 



65 ' 30375 



The computation results of M 0Q (h) for h G (/too, 00) and Mm (ft) for h G (— 00, /iio) are 
shown in Figs. |9ta) and[9]^b), respectively. As shown in Fig. [9]^a), the sign of M 00 (h) agrees 
with that of fioo < for < h — /ioo 1, and the function Mqq{K) crosses a critical value 
at h = h* 5 e (12.6197809949, 12.6197809950) at which it changes sign. Figure Mh) shows 
M w (h) for h G (— 00, h w ), whose sign agrees with that of H12 < for < h w -/i<l. 
Moreover, M 10 (h) crosses a critical value at h = h* G G (-3.1388150376, -3.1388150375) at 
which it changes sign. Therefore, for this case, in addition to the two small limit cycles, there 
also exist at least two large limit cycles bifurcating from the two different closed orbits L^* 
and Lh* of (??). One large limit cycle surrounding the center (0,0) is shown in Fig. HPT a), 
while another large limit cycle enclosing the center (1,0) with 2 small limit cycles is depicted 
in Fig.QZKb). 

It is noted that all the four sets of values of 01 and a 4 chosen above in (A), (B), (C) and 
(D) satisfy 

ai + 2a4 2n 

= — , where n is an integer and m is an odd integer, /o82 (4.5) 

ax m 

so that a consistent integrating factor (and so a consistent Hamiltonian function for the 
whole transformed system) is obtained. However, this condition is not necessary since the 
singular line 1 + a\ x = divides the phase plane into two parts, and the analysis does not 
need the continuity on the singular line. To demonstrate this, in the following we present a 
case for which the condition (??) is not satisfied. Consider the (2, 0)-distribution, and choose 
a\ = — 5 and 04 = —4. The point (ai,a 4 ) = (—5, —4) is marked by a dark circle in the 
a\-a± parameter plane (see Fig. [p. Then, 

a x + 2a 4 _ 13 26 

— ~Ti O01 — aio, On — "7T a 10> 

ai 5 3 

and 
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H(x,y) 



x 2 + y 2 1 
V h e (0, oo), when x < -, 



2(l-5x) 8 / 5 ' v ' " 5 

''" + /r V /; c (-oo, - i 2 4 / 5 ), when x > 1 



2(l-5x) 8 / 5 ' v ' 32 " 5' 

For this case, ^ 2 and /x 10 become 

M02 = ^-7rai <0 and /xi = - — tt ai < 0. 

The computation result of M 00 (h) shows that M 00 (h) < for < h <C 1, agrees 
with the sign of /io2- Moreover, Moo(O.l) = 0.0510077880 > 0, implying that there exists 
h = G (0, 0.1) such that M o(/4) — 0, and so a large limit cycle bifurcates from the closed 
orbit L h * of (??). The result of M w (h) also shows that M 10 (h) < for < - ^2 4 / 5 -/i< 
1, agreeing with the sign of (i 10 , and that M w {- ^ 2 4 / 5 -0.8) = 7.4630743072 > 0, implying 
the existence h = h* 8 e (- ^ 2 4 / 5 - 0.8, - ^ 2 4 / 5 ) = (-0.8544094102, -0.0544094102) such 
that M 10 (/ig) = 0. Thus, there exists another large limit cycle bifurcating from the closed 
orbit Lh* of (??). Therefore, this case exhibits 2 small limit cycles and 2 large limit cycles, 
leading to the existence of at least 4 limit cycles. 

Summarizing the above results with the continuity of parameters a± and 04 shows that 
at least for some regions in the ai-a 4 parameter plane the reversible near-integrable system 
(??) can exhibit at least 4 limit cycles around the two singular points (0,0) and (1,0) with 
distribution ether (3, 1) or (1,3). 

The proof of Theorem 4.1 is finished. □ 



5 Conclusion 

In this paper, we have proved that a quadratic non-Hamiltonian integrable system with two 
centers can have at least 4 limit cycles under quadratic perturbations, with distributions 
either (3,1) or (1,3). This result gives a new record, answering the open problem of the 
existence of limit cycles in near-integrable quadratic systems. It is shown that such systems 
can have at least 4 limit cycles for some regions in the 2-dimensional parameter plane, 
associated with the parameters of the integrable systems. Further research is needed on 
global analysis for all possible parameter values in the parameter plane. 
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